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Q^l Abstract 

r-| ! The Schwinger-Dyson equation for the quark self-energy is solved in the quenched laddei 

approximation for several cases of one- and two-quark-generations. The exchanges of stan 

XI 

;_! ■ dard model gluons and Higgs bosons are taken into account. It is found that Higgs bosoi 

exchange dominates the quark self-energy in the ultraviolet region for sufficiently large in 
put quark masses (> 75 GeV), causing the running quark propagator mass to increase witl 
energy-scale. The running of the quark mixing angles is also considered. No running of th( 
quark mixing angles is found for input quark masses up to and including 500 GeV. 



I. Introduction 

Unlike the renorinalization group equation, the Schwinger-Dyson equation (SDE) analysis 
presented here enables us in principle to calculate the running quark mass functions fo: 
very light quarks at low energies where quantum chromodynamics (QCD) is in the non- 
perturbative region. This is of interest because experimentally accessible quantities an 
presently found at low energy-scales. In the presence of heavy quarks, such as bottom anc 
top, it allows the calculation of running mass matrix at momentum scale of the order o: 
less than the quark masses, at which scale the RG equation results may not be reliable 
Moreover, this nonperturbative SDE approach is also interesting in that it can describe th( 
running mass function of an ultra-heavy quark, with a large Yukawa coupling that canno' 
be treated reliably by perturbation theory. This would be relevant for a fourth generatioi 
quark, for example. 

The calculation of the variation of mixing angles with momentum scale is of great ex- 
perimental interest since they are usually measurable only in a very limited kinematic re 
gion. To illustrate this point, consider the top-strange CKM mixing parameter, Vts- Direc 
top-antitop quark bound state production at the next linear collider would occur at approx- 
imately > 300 GeV^, and Vts would be evaluated there. Presently, some estimations of Vt, 
are being done indirectly using penguin diagrams in the calculation of 6 — > s transitions 
In this method a virtual top quark mixes with a virtual strange quark at a few GeV^ ||I[ 
Clearly if Vjs(300 GeV^) 7^ Vts{— 3 GeV^), complications would arise. This possibility 
should be investigated both by RG techniques and by the SDE. In particular SDE studie; 
can reveal the existence of any nonperturbative effects which might be absent in RG studies 

There have been several studies of the quark propagator Schwinger-Dyson Equation 
These have fallen into the areas of chiral-symmetry breaking p|, |^, ^ |^, ^, confinemeni 
[[^], and related phenomenology including pseudoscalar decay constants and form factori 



[P, |y, |10|, [TT| . The complexity of the SDEs makes it neccesary to use various approximation; 
or assumptions. The numerous QCD quark SDE studies make assumptions about the quark 
gluon vertex and they model the effective gluon propagator. Many of these studies utiliz( 
the Landau gauge, ladder approximation and the angle approximation to express the Integra 



SDE in terms of a differential equation [0 . Recent reviews of SDE literature which include 



discussion of and references for the above points include [|12], |T3], . 

Though the propagator is a gauge dependent object, and this is true in the approximation! 
that we use, our results do agree with the gauge-invariant one-loop renormalization grouj 
analysis results in the ultraviolet region. We are careful not to neglect derivatives of th( 
strong coupling, asiq"^), to be consistent with the one-loop RGE result [p!5[ . 



Roberts and McKellar |T^ have studied the validity of the angle approximation used ii 
the context of the model of Atkinson and Johnson |^. They use the Landau gauge anc 
the running strong coupling ^^(g^) = dn/ln^xo + q'^/Aq^j-,), and re-express the Integra 
equation as a differential equation. Roberts and McKellar find that the angle approximatioi 
is qualitatively reasonable, introducing errors of order 12%. For other models, however 
the angle approximation is not so reliable, particularly when the strong coupling become; 
singular at the origin. To minimize problems with the angle approximation, we utilize som( 
aspects of the model of Atkinson and Johnson p. 

Rather than investigate chiral symmetry breaking, this study assumes that chiral sym- 
metry is broken and solves the quark propagator SDE to find running quark mass function; 
and running quark mixing angles. In addition to the usual gluon contributions, we include 
Higgs boson exchange effects in the quark self-energy. The contributions to the quark self 
energy due to Higgs bosons are being studied here for the first time in the context of th( 
quark propagator SDE. We also study for the first time in a Higgs-boson-plus-gluon mode 
the multigenerational cases of quarks which enable us to analyze quark mixing angles. 

II. Model 

The general form of the quark propagator SDE that we use is 

S-\q) = S,\q) - / -0^,9sl,Sik)gA.G^-^ - / ^,gHS{k)gyAj,PH, (l 

where (7^7^ and gsAi, are quark-gluon vertices, S{k) is the quark propagator, G^'^ is th< 
gluon propagator, gyAn is the quark-Higgs boson vertex factor, and Ph is the Higgs bosoi 
propagator, to be discussed below. 



In this study we will simplify the quark propagator by using the ladder approximation 
namely we will set all the vertices to be the bare vertices. We will, however, use the RG 
improved expressions for both the QCD and Yukawa couplings. In the Landau gauge A{q) = 
1 in the approximate form used here of the Schwinger-Dyson-Equation (SDE) in QCD ||T7| 
Thus the ladder approximation quark propagator in the Landau gauge in Minkowski spac< 
that we use is iS{k) = {!/+ M){k'^ — M^)~^, where M is a 2 x 2 quark mass function matrix 
assumed to be symmetric for this study. The ladder approximation quark propagator if 
accurate enough for our purposes because this is an exploratory study of a quark self-energ] 
SDE which includes QCD and Yukawa contributions [|T8| . 

The gluon propagator we use is 

-9'" + -^j G{k) (2 

with 

and the leading-log (one-loop) coupling 

9sik) = — ^T- (4. 

ln{xo - j^) 

where d = 12/(33 — 2nf), Uf is the number of flavors. This form has the correct — /c^ -^ oc 
leading order QCD behavior (for — A;^ > 0). This is a simple extension of the form usee 
successfully to fit J/ip and T spectroscopic data [|T^. This leading- log QCD coupling if 
renormalization-scheme and gauge independent, as is its second order, two-loop, extensioi 
[ pO[] . xq is a parameter that functions as a smooth infrared cutoff. Equation (^) does no' 
accurately model the gluon potential in the infrared region, but it does enable us to assess 
the relative magnitude and shape effects of the Higgs boson and gluon exchanges in th( 
quark self-energy [pl|] . 

The QCD coupling also contains the parameter Uf, the number of quark flavors. Th( 

2 

number of quark flavors is determined by the energy scale, y = -t4 — , where q is the momen 
tum transfer. We use a step-function expression to increment Uf at q =L27 GeV (n/ = 4) 
4.25 GeV (n/ = 5), 160 GeV (n/ = 6), and 500 GeV (n/ = 8), where the running quarl 



mass values mc{m^ = 1.27 GeV, and mb{mi,) = 4.25 GeV are from Gasser and Leutwylei 



22 



We must also model the Yukawa couplings. We consider a 2 x 2 symmetric Yukawa 
coupling matrix, which can be diagonalized by a 2 x 2 unitary matrix, to give diagona 
elements g±. The one- loop Renormalization Group Equation for the diagonalized Yukawc 
couplings, g±, without weak interaction effects, is 

dg± 1 /9 



9iit)-^9QCDit)9±it)) (s; 



dt 47r2 V8 
with t = ln{xo + q^ / IS^QCj;)) /2. This has solution 

9±- ^ .2d 9 t_ 1° 

^^'^ 167r2 l-2d 

where the C± must be determined by boundary conditions. We apply the boundary con 
ditions m'^{l GeV) = g'^{l GeV)v^/2 for small (< 150 GeV) running quark masses, anc 
m^(m-i-) = 5'±(^7i±)v^/2 for large running quark masses, where v ~ 246 GeV, and m+ = m 
and m_ = rric, if we consider the 2x2 mass matrix to represent the second and thirc 
generations, for example. We chose Aqcd = 0.18 GeV to set the boundary values of g±. 

We convert Eq. (|1]) into 4-d spherical polar coordinates, and integrate over and } 
analytically. To integrate over 6, we employ the so-called angle approximation |T5[ in th( 
QCD sector and its analog in the Yukawa sector. The angle approximation is equivalem 



to keeping the first term in an expansion of the angular, 6, integral |2i|. It predicts th( 
same leading behavior for the ultraviolet asymptotic forms of M(y) as given by operator 
product-expansion analysis in the pure QCD case |25|. Roberts and McKellar [|1^ showec 
that the angle approximation may be useful for a qualitative study of the quark SDE whei 
the infrared behavior is expected to be smooth, as in our gluon model. We also make th( 
approximation Mjj -C Aq(^,jj{x + y — I^Jxycosff). This approximation is appropriate for th( 
heavy quark cases that we study. For the light quark cases, this approximation exaggerate; 
Higgs effects, thereby providing an upper bound to any new nonperturbative effects. For al 
input quark masses this Yukawa sector approximation holds in the asymptotic region wher( 
y > 10^ 



The resulting 1-d SDE is thus 
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Mq + AqcdJ xdx 
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We solve Eq. (|^ by converting it into a second order differential equation (see Appendix 

and applying a fourth-order Runge-Kutta subroutine. To ensure that the differential equa- 
tion is equivalent to our integral equation, we must enforce some initial conditions. Th< 
QCD quark propagator integral SDE relates the values for the derivatives of the quarl 
masses at the origin to the quark masses themselves. For the pure QCD 1-particle case 
"'^^"° , where aQcoi^) = 'iTd/ln{xo) [^. The initial conditions for our QCE 



dy \y=o 



2A/(0)7i 



plus Higgs boson study must be consistent with this. For the two-quark-generation case 
we assume the propagator quark mass function takes the form Mij{y) = Mqij + Mnj y nea: 
the origin, where M^ij are the initial input quark masses, Mg(0), and the forms of the firsi 
derivatives of the mass matrix elements, Mnj, are given in the Appendix. The pure QCE 
sector tells us further that since we have included a nonzero bare quark mass (see Eq. ^, w< 
must agree with the irregular QCD mass solution in the asymptotic region, i.e. for y ^ 1 
M{y) oc {ln{y + Xq))~'^ as long as QCD is the dominant effect [|5|. We will use this as ai 
asymptotic test of our numerical results. 



quark 


M,(0 GeV) 


mg{l GeV) 


up 


560 MeV 


5.6 MeV 


down 


560 MeV 


9.9 MeV 


strange 


720 MeV 


199 MeV 


charm 


1.66 GeV 


1.32 GeV 


bottom 


4.89 GeV 


4.52 GeV 



Table la: Input quark masses 



Numerical input values are given in Tables 1. 



quark 


M,{0 GeV) 


mq{mq) 


top 


179 GeV 


17 A GeV 


bottom' 


506 GeV 


500 GeV 


top' 


506 GeV 


500 GeV 



Table lb: Input quark masses 

We use two sets of input quark masses. For up, down, charm, strange, and bottom quark; 

values as initial conditions on the Mg(0) for th( 



at q=0 we use Jain and Munczek's 



integration subroutine. The 1 GeV values are used to fix the values of C± of the Yukawa 
sector (except for the top quark and fourth generation quarks). The light quark (u,d,s) '. 
GeV input masses are obtained by choosing values at 1 GeV which are consistent with th( 
1994 Particle Data Book and Y. Koide's paper "Table of Running Quark Masses" [^. Th( 
heavy quark (c and b) 1 GeV mass values are obtained by running the heavy quark mas! 
values given in the 1994 Particle Data Book to 1 GeV |^|. The Schwinger-Dyson evolutioi 
of M(g) is small below Tng{q) for heavy quarks |^, as seen from the slightly higher evolvec 
values of Mt^t',b'{0) starting from mq{mq) values. For the top quark and the fourth generatioi 
quarks we normalize C± at m±{m±) as indicated after Eq. @. 

For the multiple-quark-generation case we must also input the g^/Ag^^ = values o 
the mixing angles for the up-sector and down-sector. We somewhat arbitrarily choose th( 
up-sector mixing angle to be 9up = 0.5 radians. To ensure that the cabibbo mixing angl* 
agrees with the accepted SM value for the first two generations, Ocabibbo = 0.22 radians a' 



small g^ [EH], we choose 6', 



down 



Oup — cabibbo- For the third-fourth two-generation case, ou: 



convention is to adopt the same Ocabibbo input value, since it is unknown. For the second-thirc 
two-generation case, we utilize the leading angle expansion term of the 1994 Particle Dat< 
Book's Vcb, sm^23 = 0.04 |^. We are primarily interested in the question of the running o 
the mixing angles, not their specific values, in any case. Next we turn to the results of this 
analysis. 



III. Results 

We solve the one-quark-generation version of Eq. (0) for seven cases of quark input masse! 
given in Table 1 for the running quark mass functions. We summarize some characteristic; 
of these results, including the initial and final values of the quark mass functions in Tabl< 
2, where A and B are the initial and final integration points, for example 5 x 10~^, anc 
9.77 X 10^, respectively, in units of A^^^. We also compare the QCD plus Higgs bosoi 
exchange case with the QCD-only case in the asymptotic region to get the ratio ?f^^" ■ 



QCD 



quark 


MM) 


M,{B) 


^Iqcd+h 

Mqcd 


up 


s.hAqcd 


I.OOAqcd 


1 


down 


s.hAqcd 


I.OOAqcd 


1 


strange 


S.99Aqcd 


IA2Aqcd 


1 


charm 


9.22Aqcd 


4.21Aqcd 


1 


bottom 


27.2Aqcd 


15.3Aqcd 


~1 


top 


978Aqcd 


825Aqcd 


1.01 


bottom'=top' 


2790Aqcd 


3140Aqcd 


7^1, varies 



Table 2: One-quark-generation results, for A = 5 x 10 , B = 9.77 x 10 



We solve the second order differential equation version of Eq. (0), with initial condition; 
given in an Appendix, to get the quark mass function weak eigenstates Mn, M12, M22. We ge' 
the quark mass function eigenstates from M± = 



Mil + M22 ± J{M: 



'11 



M2^^2 



'22 



4M2 



12 



and the mixing angles from tan2ef = 2Mu[{M22 - Mn) ± J{M22 - Muf + AM\ 



li-^ ^^ when 



/ stands for the flavor-sector, and where OcaUhho = Gup — (^down- We look at three cases o 
two-quark-generations: first and second generations, second and third generations, and thirc 
and fourth generations. 

In Figure 1, we graph the bottom quark mass function for the QCD plus Higgs boson cas< 
(solid line) and for the QCD-only case (dashed line) from the second and third generation; 
analysis. The two lines are coincident. We also graph the top quark mass function for ou: 
QCD plus Higgs boson model (solid line) from the second and third generations analysis ii 
Figure 2, with the corresponding QCD-only case (dashed line) shown for comparison. Th< 



QCD plus Higgs boson case of the top quark mass function shows some deviation from pur* 
QCD. The Higgs boson term has an opposite sign relative to the gluon term, so it drives th< 
mass function up, in contrast to the usual QCD-only decreasing mass function. In Figure I 
we illustrate the dependence of the fourth generation top' quark mass function on the energy 
scale, from the third and fourth generations case. The Higgs boson has a substantial effeci 
on the fourth generation quark self-energy, i.e. the QCD plus Higgs boson case (solid line 
disagrees substantially with the pure QCD case (dashed line). In fact, the running top' quarl 
mass is increasing in the asymptotic region. In Figure 4 we graph the 'running' Cabibbc 
mixing angle versus the energy scale for the third-fourth generations case, which shows th< 
largest Higgs effect on the mass. There does not appear to be any variation in the mixing 
angle with respect to energy scale. This is true in all cases we study. 

In Table 3 we summarize the two-quark-generation results. The two-quark- generatioi 
results are very similar to the one-quark-generation results, as we see by comparing Table; 
2 and 3. Apparently the mixing effects in our model in the two-generation cases do no' 
significantly affect the quark mass functions. 



quark 


MM) 


M,{B) 


Mqcd+h 
Mqcd 


up 


s.hAqcd 


I.OOAqcd 




down 


s.hAqcd 


I.OOAqcd 




strange 


S.99Aqcd 


IA2Aqcd 




charm 


9.22Aqcd 


^■21Aqcd 




bottom 


27.2Aqcd 


15.3Aqc'd 




top 


993Aqcd 


840Aqcd 


1.01 


bottom'=top' 


2810Aqcd 


^IbQAQCD 


7^1, varies 



Table 3: Two-quark-generation mass results, for A = 5 x 10 ^, B = 9.77 x 10^. 



IV. Conclusions 



All calculational results have been verified by an independent calculation as indicatec 
in the Appendix. We have calculated the momentum-dependent, quark propagator mass 



for seven flavors of one-quark-generation quarks via a nonperturbative SDE treatment witl 
Higgs boson exchange contributions. We do not see the effect in the ultraviolet region on th( 
quark self-energy due to adding the Higgs boson interaction until the input quark mass if 
sufficiently large, Mg{0 GeV) > 75 GeV. Thus only the top quark and the fourth generatioi 
quark results differ from pure QCD results. 

We have also calculated the quark mass function for eight flavors of two-quark-generatioi 
quarks via a nonperturbative Schwinger-Dyson matrix equation analysis with Higgs bosoi 
exchange contributions to the quark self-energy. The results are similar to the one-quark 
generation results. Only the top quark and fourth generation quarks, top', and bottom' 
differ signiflcantly from the pure QCD results, with the fourth generation quarks showing 
the most marked effect of the Higgs. For the fourth generation quarks, the quark mass 
function actually increases with energy-scale, due to Higgs boson exchange. 

We can analytically determine where the Yukawa term dominates. In terms oi y = 
q^/AqcD and d = 12/(33 - 2n/) we get 

9d \ (2)2^ 



y> 



exp 



Xo 



2d- I J 32TT^dC±_ 
where C± < ioi y > —xq + exp{ ^*^'^/ ^^J^}- Thus for the top quark, Ct = —0.19, anc 



the Yukawa term dominates when q is beyond the Planck scale. For the fourth generatioi 



quark, C^th = —0.053, and the Yukawa term dominates when q > 1.5 TeV. 

We have calculated the running mixing angle, which we call theta-cabibbo, for thre( 
cases of two-quark-generations. We detect no running of theta-cabibbo or any other mixinj 
angle due to any effect, for any quark mass input value up to and including 500 GeV. Thif 
is expected for pure QCD and can be proven analytically [^ . This provides a check of th( 
numerical work. By studying the mixing angle results, we determine that the running of th( 
mixing angles due to nonperturbative effects in our model is less than one part in 10^, far less 
than any conceivable experimental sensitivity. We have been unable to prove analyticalh 
that the cabibbo angle is constant when the Yukawa interaction is included in the SDI 
analysis, but we have established numerically that there can be no observable effect ever 
for the heavy hypothetical fourth generation. Based on our results, it appears to be safe t( 
ignore running of Vts and Vtb, for example, in extrapolating from ml to mf in parametrizing 



data [|1|, ^. Our nonperturbative results agree in this respect with the RGE analysis of th( 
model that we study, since the RGE's can always be diagonalized completely and no running 
of the mixing angles occurs. As remarked earlier, the ^/-dependence of our mass function; 
agrees with that of the 1-loop RGE's in the asymptotic region, though not in the region a 
or below the scale of the input masses. 

In our investigation we have neglected the contributions of Goldstone bosons, W^, 1? 
and 7 exchange in the quark self-energy. Of these, we expect only the Goldstone bosons t( 
have a large effect on the quark self-energy in the Landau gauge, and then only in the heavj 
quark sector where it is strongly coupled. Thus the next step in this analysis is to expanc 
our model to include Goldstone boson exchange. For heavy quarks, we expect the effects or 
the quark self-energy from the electroweak bosons to be significantly smaller than those o 
the gluon, Higgs boson, and Goldstone bosons. 

We could also expand the flavor sector to include the three-quark- generation case anc 
possibly the four-quark-generation case. This would enable us to calculate the running o 
the elements of the three-quark-generation Cabibbo-Kobayashi-Maskawa matrix and its four 
quark-generation analog. These calculations should be undertaken with due consideratioi 
however, since introducing additional quark generations significantly increases the length o 
the calculations. Moreover, our two-quark-generation case seemed to indicate that the effect; 
on running quark mass functions due to mixing between generations are negligible. 

This method is well-suited to study a variety of nonperturbative non-Standard-Mode 
effects. For example, we could investigate the case of the third generation top quark ex- 
periencing a new gauge interaction instead of, or in addition to, the Higgs mechanism b] 
replacing the Higgs boson exchange with a new gauge boson exchange in the top quarl 
self-energy. This would enable us to study models of dynamical breaking of electroweal 
symmetry where the top quark plays a special role [^, |3^ . Alternately, we could study th( 
case of the fourth generation quarks experiencing a new gauge interaction, in addition to th< 
Higgs mechanism, which the SM quarks do not experience, to study new, dynamical effects 
We address these and other, related, points in a future publication. 
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Appendix 

Using g± we can specify giji 

gl, = g^s^n-e + g^cos^O, (S; 

fi'12 = {.9+ - 9-)cos9 sine, (9 

(^22 = g+cos^9 + g^sin^O. (10 



Then in terms of 



and 



A = ^ ^" ^ ^^^ (11 

yHn{xQ + y) 32n'^y'^ 

„ _ 912(911 + 922) ,-.^ 

^ ~ 327r2t/2 ' ^^^ 

C = ^ ^12 + 9I2 ^3 

yHn{xQ + y) 32-K^y^ 



V = {yAlcD + Ml^ + Ml^){yAlcD + Ml^ + Ml^) 

- MUMu + M22f, (14 

Mn = Mn{-yAlcD-M^2)-M^2M22, {15 

Mu = MnMuM22 + Mu{-yAlcD~Mf^), (16; 

M22 = Ml^Mu + M22{-yAlcD-Ml^), (17; 



we can write the components of the second order differential equation: 

£Mii _ AqcdV, 



where 



dy' 



V 



-{MuA-MuB) 



+ 



1 



^^ {SC + TB) + --^[TC + QB) 



AC-B'^\ dy 

(PMi2 _ AqcDy 



dy^ 



V 



dy 
-MuB + MuC) 



+ 



1 






d'^M22 ^ Aqcpy 
dy^ V 



-M12B + M22C) 



+ 



AC-B^ \ dy 



[SB + TA) + —-^[TB + QA) 



dy 



£ 



2A d 

y yHn'^{xQ + y){xQ + y) 



2S 
T = — 

y 



d 



y yHn?{xQ + y){xQ + y)' 
The initial condition Mnj are given by: 

dMn(|/), _ ^QCDy'' 



dy 

dMujy) 

dy 
dM22{y) 



\y=0 



\y=o 



dy 



\y=o 



2-V 

^QCDy\ 

2-V 

^QCDy\ 

2-V 



{MiiA-Mi2B) = M^^^ 



-MuB + M12C) = Miu 



-M12B + M22C) =Mi22- 



(1^ 



In [^ we use S ^ = —i[(/A — B] and solve the integro-differential system 



A-\q') = l + ^J^^k-qgy[e-M\k^)rgy 



1 



q'J (2 






:i + 0A;-g + 2(e-i; 



{k - qf - Mfj 

q- {k — q)k - {k — q) 
{k-qf _ 



.[e-M\e)r^^^ (27: 



M{q')=M,A-^ + ij^^gYM[e-M\e)]-'gy 



(27r)4 ^" L "^ " ^' {k - qf - Ml 



where ^ is the gauge parameter. Our gauge choice here is ^ = 0. Equations ( p7D and (pBD ar* 
gauge independent at large momenta where the terms depending on ^ are subleading. W< 
have solved Eqs. (|27|) and (pS]) by using the angle approximation for the QCD terms, taking 
gy as a constant, and performing the angular integrals exactly for the Higgs term with i 
massive Higgs boson. The result is converted into integro-differential equations and solvec 
numerically. The solutions for the mass functions agree very well with solutions to Eqs 
(p!8D-(pOD for a variety of boundary conditions. Further applications of the A ^ 1 solution; 



are presented in ||29 
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Figure Captions 

Figure 1: mhottom{y) / ^qcd vs. y for A = 5 x 10^^, B = 9.77 x 10®, and inputs xq = 5 
mo =4.89, 506 GeV. Tlie solid (daslied) line is the gluon plus Higgs boson (QCD-only 
interaction result. 

Figure 2: rritopiy) / J^qc d vs. y for A = 5 x 10"^, B = 9.77 x 10®, and inputs xq = 5 
rriQ =179, 506 GeV. The solid (dashed) line is the gluon plus Higgs boson (QCD-only 
interaction result. 

Figure 3: mtop'{y)/ J^qcd vs. y for A = 5 x 10"^, B = 9.77 x 10®, and inputs xq = 5 
mo =179, 506 GeV. The sohd (dashed) line is the gluon plus Higgs boson (QCD-only 
interaction result. 

Figure 4: OcaUbboiy) vs. yioT A = 5x 10"^, B = 9.77x 10®, and inputs Xq = 5, mo =4.89, 506 
179, 506 GeV. The solid (dashed) line is the gluon plus Higgs boson (QCD-only) interactioi 
result. 
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